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> ; Abstract 
oo 

' We derive an approximate relic neutrino asymmetry evolution equation that 

■ systematically incorporates repopulation processes fi'om the full quantum ki- 

. netic equations (QKEs). It is shown that in the collision dominant epoch, the 
said equation reduces precisely to the expression obtained previously from the 

CL|| static/adiabatic approximation. The present treatment thus provides a rig- 

Q^' orous justification for the seemingly incongruous assumptions of a negligible 

^ ■ repopulation function and instantaneous repopulation sometimes employed in 

• • . earlier works. 
• 1— I - 

H ■ I. INTRODUCTION 

The computationally convenient and rewarding static approximation [0 may be rigor- 
ously derived from the quantum kinetic equations (QKEs) [Q,^ in the adiabatic Boltzmann 
limit 0,^, assuming that the repopulation or refilling function is insignificantly small. In ac- 
tual numerical studies, however, we also adopt the view that a depleted neutrino momentum 
state is instantaneously refilled from the background plasma. The latter condition naturally 
implies that while the repopulation function may be small, it cannot be identically zero. 

In principle, simultaneous suppositions of a vanishing refilling function and instantaneous 
repopulation are somewhat contradictory, or at least uncontrolled. However, their combined 
validity in the collision dominant epoch is strongly supported by numerical evidence [|T],0. 
Indeed, instantaneous repopulation without the assumption of a negligible repopulation func- 
tion can be consistently implemented as a sensible approximation to the QKEs 0. 

Commencing with the full-fledged QKEs, we demonstrate rigorously that the two afore- 
mentioned approximations, conceived originally for computational convenience, are indeed 
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appropriate in a certain limit. Our approach differs from the fully consistent instantaneous 
repopulation approximation introduced in Ref. since we do not make any a priori as- 
sumptions regarding the repopulation rate, instantaneous or otherwise. Instead, the present 
work shows that a "brute force" treatment of the QKEs incorporating a finite repopulation 
rate does indeed lead to approximate evolution equations that reduce to those following from 
the dual approximations of a negligible refilling function and instantaneous repopulation in 
the temperature regime of interest. 



II. NOMENCLATURE 

Consider a two- flavour system comprising an active species z/q, (where a = e, yU or r), 
and a sterile species z/^, whose properties are characterised by the density matrix |^,^ 

p(p) = ^[Po(p) + P(p)-a], (1) 

in which the variables Pq and P(p) = Px{p)x + Py{p)y + Pz{p)z are functions of both time t 
and momentum p, and a = axX + ayij + OzZ are the Pauli matrices. The diagonal entries of 
p represent respectively the and Vs distribution functions in momentum space, that is, 

Ns{p) = llPoip) - nb)]iV^^(p, 0), (2) 
where the reference distribution N'^^{p, 0) is defined as the equilibrium Fermi-Dirac function, 

with chemical potential /i = at temperature T. 

The evolution of Po{p) and P{p) is governed by the quantum kinetic equations (QKEs) 

(9P BP 

— = V{p) X P{p) - D{p)[Pxip)x + Py{p)y] + -^z, 

riP 

-gf = RM- (4) 

Here, the quantity V(p) = (3{p)x + \{p)z is the matter potential vector J^, with 

/3(p) = ^^sin2^o, 
2p 

Kp) = -^IKp) - <P) - cos2^o], (5) 

in which Am^ is the mass-squared difference between the neutrino states, 9q is the vacuum 
mixing angle, and 
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, , 4C(3)v^G'^LHT3p 
a P = ^T-^ , 

4C(3)V2GfA.TV 
^(^^ = vr^Am^m^^ ' 

given that is the Fermi constant, my/ the VT-boson mass, C the Riemann zeta function 
and Af. ~ 17, v4^^^ ~ 4.9. The effective total lepton number (for the a-neutrino species) 

L(") = L^ + L, + L^ + Lr + r] = 2L^ + L, (7) 

combines all asymmetries individually defined as = {fii,^ — n-p^)/n^, where the symbol n 
denotes number density, and rj is a. small term due to the cosmological baryon asymmetry. 
The decoherence function D{p) is related to the total collision rate for z/q, T{p), via |^ 

D{p) = Inp) = \-^y^GlT\ (8) 

where (p)o — 3.15T is the average momentum for a relativistic Fermi-Dirac distribution 
with zero chemical potential, ye — 4, y^^,- ~ 2.9 and Ze — 0.1, z^^r — 0.04. 
The repopulation or refilling function 

Ra{p) ^ T{p) - i[Po(p) + PM]] , (9) 

with 

is determined from assuming thermal equilibrium for all species in the background plasma, 
and that the z/q, distribution is also approximately thermal 0]. Physically, Eq. @ means that 
as active neutrinos of some momentum are converted to sterile neutrinos, a gap is created 
in the Fermi-Dirac distribution. The weak interaction processes repopulate this depleted 
state, driving the ensemble back towards equilibrium. 

A separate but equivalent set of expressions parameterises the evolution of the Vs 
system. These are distinguished from their ordinary counterparts with an overhead bar. 



III. STANDARD STATIC/ADIABATIC LIMIT: A BRIEF OUTLINE 

Together with the requirement of a + s lepton number conservation, we may derive from 
the QKEs an exact evolution equation for the neutrino asymmetry Lq, that is 

^ = ^ / m{p) - Pyip)Wip, mp. (11) 

where the quantities Py and Py are obtained from numerically integrating the respec- 
tive QKEs for the neutrino and antineutrino systems given by Eq. (^. The role of the 
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static/adiabatic limit approximation therefore, is to generate approximate expres- 

sions for Py that are dynamically decoupled from the other variables in order to minimise 
the computational effort. 

The first step in the formal adiabatic procedure of Ref. consists of setting the repop- 
ulation function to zero, i.e., 

Ra ^ 0, (12) 

thereby reducing the four-component QKEs to a system of three coupled homogeneous 
differential equations. Further manipulation produces the approximate equality 



D2 + A 

in the limit D, \X\ ^ such that Eq. (|1T]) become^ 



dLa 1 
dt 2n^ 



dp. (14) 



Detailed discussions of the standard adiabatic procedure may be found in the companion 
paper Ref. p. 

In numerical studies, an analogous expression describing sterile neutrino production for 
each momentum state is employed for tracking the quantity Ns{p) in Eq. (p!^, while the 
distribution function is taken to be 

iV„(p)~iV'=i(p,^). (15) 

This is the so-called instantaneous repopulation approximation, which assumes that a de- 
pleted momentum state is immediately refilled from the background medium so that thermal 
equilibrium is always maintained. 

A naive interpretation of Eqs. and (0) suggests that instantaneous repopulation 
leads directly to an identically zero refilling rate, which seems to argue for consistency 
between the central assumptions contained in Eqs. (|l^) and (|lD. However, the instantaneous 
repopulation limit is also related to taking the collision rate to infinity. Thus the right hand 
side of Eq. (^) is an a priori undetermined finite and generally nonzero function [Q, that is 
apparently at odds with the assumption of a vanishing repopulation rate. 

IV. EXTENDED ADIABATIC APPROXIMATION 

The exact QKEs in Eq. (^ are fundamentally a system of four coupled first order differen- 
tial equations, of both the homogeneous and inhomogeneous varieties, that is best displayed 
in matrix form. 



^The denominator in Eq. (|T^), and therefore Eq. (14), should properly contain an additive term 
P'^ to ensure the equations' validity even when the condition D ^ does not hold (see companion 
paper Ref. |5|). In this work, however, we shall assume that the said condition is always met, and 
omit the /J^ term for simplicity. 
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di 



(P.\ ( 



Pz 

\PoJ 



-D -X \ 

X -D -13 

p -D -D 

V -D -D J 



(P.\ ( \ 



P^ 
P. 

\PoJ 





2DK, 



= IC'^P'P + A. 



(16) 



Dealing firstly with the homogeneous part of Eq. ([I6|), we introduce an instantaneous diag- 
onal basis onto which we map the vector P^'p = (P, Pq) from its original fixed basis via the 
transformation 

/ \ 

= Q^P = W.pP'-P, 



\Q7J 



(17) 



where, by definition, 



ICJ = diag (Ai, A2, A3, A4) = Ur^U^U, 



-1 



(18) 



such that the matrix K,']' is diagonal and similar to /C^p. 

The eigenvalues Aj are solutions to the quartic characteristic equation 

A^ + ADA^ + (5L>2 + + f3^)A'^ + 2D{D'^ + + p'^)K + = 0, 



(19) 



formally given by 

^Jd^ - A2 - /32 ± + 2D^A^ + A^ + 2/3^A^ - 2p'^D'^ + 
Aj — —D ± — ^= . (20) 

We are primarily interested in the limit D, \X\ ^ \f3\, in which case the eigenvalues are well 
approximated by the leading order terms in the small f3 power series expansion of Eq. (pUD, 

P^X 



A =JJ^ 
^ 2 + A2 



+ 0(/3' 

?2 



A4 = -2D + ^^^f^ + 0(/5^). (21) 

column eigenvectors which for A 7^ are 



The transformation matrix U^^ comprises the 
exactly 

/ 1 \ 

\Ai{2D+Ai) 
l3D{D+k,) 
\ AA,(2D+A,) / 

while its inverse Wrp consists of the row vectors 



rp 



(22) 
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rp_ 



X 



n 



PD A, - A, 



AD^ + 2D\Ki + A„ + A„) + D{KiK 

+ 2D(A; + A„ + A„) + AiA^ + A^A„ + A„A, 



A/AmA„ 



(23) 



where the indices /, m and n are the three integers not equal to e.g., for i = 1, I, m and 
n are 2, 3 and 4 respectively. 

Note that Eqs. (^) and (^31) do not apply if two or more eigenvalues are degenerate, in 
which case the matrix /C^ may genuinely have less than four distinct eigenvectors thereby 
rendering momentarily uninvertible, or, if four linearly independent eigenvectors exist, 
alternative methods are needed for their evaluation. However, such circumstances are hard 
to come by and have practically no influence on the outcome of this work. Thus the exact 
QKEs [Eq. (|16])] are equivalently 



dQ 



rp 



dt 



(24) 



written in the new instantaneous basis where the transformation matrices are valid. 

The term U^^^^ in Eq. ( ^4|) contains explicit dependence on the derivatives of the 
parameters A and (3. These are discarded in the adiabatic limit (see Ref. for the rele- 
vant constraint^) so that the remaining terms form a set of four decoupled inhomogeneous 
differential equations given by 



rp 



dt 



VfQ^^ + B, 



where B = WrpA, or equivalently in index form, 

^QrW^A,(t)gjp(t) + i?,(t), t 



1,2,. ..,4. 



Equation ( |26D may be formally solved to give 

g^P(t) = e/o ^'(*')'^*'QP(0) + e/'^»(*')'^*'j,, 
with the inhomogeneous segment of the QKEs contained entirely in the integral 

J, = f\-r'^^^'"^'''"Bi{t')dt', 

Jo 



(25) 



(26) 



(27) 



(28) 



^The bounds in Ref. |§] are calculated assuming ~ 0, and should naturally be different from 
those one would obtain with the incorporation of a finite repopulation function. However, we 
expect the difference to be minimal given the similarity between the relevant entries in the two 
cases' respective transformation matrices in the appropriate hmit. 
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thus amounting to the following solution for the vector P''^, 

/e/o^^{*')'^*' 

e/o^^(*')'^*' 

e/o^^(*')'^*' 



Py{t) 
V Po{t) J 



u-\t) 



V 



+ 



e/o^*(*')'^*' 

fePMindt' ^ 

e/*^^^*')'^*' 
e/'^^(*')'^*' 

e/*^^(*')'^*' 



vx4y 



fQ?io)\ 
Q?{0) 

\Q7{0)J 



(29) 



in the adiabatic limit. 

Observe in Eq. (^) that the real components of the eigenvalues Ai, A2 and A4 are of the 
order D, while A3 is proportional to This implies that the exponentials exp[/Q Ai{t')dt'], 
where i = 1,2,4, are rapidly damped relative to the "decay" time scale of their A3 coun- 
terpart in the homogeneous part of Eq. (|^) (for a full discussion, see Ref. P). We may 
therefore implement in Eq. (123) the collision dominance approximation 



(30) 



to obtain 

Py(t) 
Pzit) 
Po{t) 



D + A 



-XA,{2D + A,f° ^°^-^AA,(2D + a/ 

AA3(2D + A3)' ^"^+^AA,(2D + a/' 



(31) 



and consequently, from combining the above expressions for Py and 



PiD + A^ 



( A,{2D + A,) A,(2D + A,) \ P{D + A,)^ J^K.jfWj (on) 
.k^A V + A3) P{D + A,) ; AA,(A, + 2D) ^ ^ 



Note that the collision dominance approximation in Eq. (|30|) is not immediately applicable 
to the indefinite integrals associated with the inhomogeneous part of Eq. (p9D since the term 
Xj contains the factor exp[— /* Aj(t')(it']. 

With the quantity Bi in Eq. (E3) evaluating to 



I3D 



.,■ Ai- Ai 

3+1- * ^ 



(33) 
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courtesy of Eq. ( p3D , Eq. (p2D simplifies to 



Py{t) 



j\D+iX)dt 



t rt' . 



D + i\ 



' ^f'(D+iX)dt" PPKq 

Jo D-iX 



(34) 



as a power series in the small expansion parameter (3. Clearly, the first term in the above 
expression is simply the standard adiabatic result given by Eq. (|l^). 

The second item in Eq. (Q) involving the inhomogeneous parameter Ka may be partially 
solved using integration by parts, 

-j\D-iX)dt' ^ f'{D-iX)dt" I^DKg 

D + iX 



f3DKq 



+ A2 



t'=t 



^- J^{D~iX)dt' PPKq 

+ A2 



t'=0 



.-f{D-tX)dt' p p\D-tX)dt 



I3D dKg d f pD 



D2 + A2 dt' 



'dt' \D^ + X^ 



dt'. 



(35) 



A sufficiently large damping parameter D again allows the approximation of exp[— /q(Z) — 
iX)dt'] 0, thereby promptly obhterating the second term on the right hand side of the 
last equality. The quantity ^ d^+x"^ ) integral turns out to be one of the elements 

constituting the 4x4 matrix W^p^^ • Consistency with the adiabatic approximation then 



calls for the setting of ^ (75!^) — such that solutions to Prp have no explicit dependence 
on the time derivatives of the parameters D, A and (3. 

Further exploitation of the powerful technique of integration by parts on Eq. (p5D in 
conjunction with the aforementioned approximations thus yields 



-f{D-iX)dt' /•* j''{D-iX)dt" f^^^^ ^i' 

Jo D + iX 

f3DKg 



Z}2 + A2 



+ 



t'=t 



N 

E 



\\D-iXj D2 + A2c/f(") 



-I t'=t 



-p-riD-iX)dt' f\f'{D~iX)dt" 



1 



N 



pD 



D-lXj Z}2 + ^2^^/(7V+l)^"^^'- 

Above the neutrino decoupling temperature, the function is well approximated by 

oP/T 



(36) 



7r2 1 + 



(37) 



for small L^,, and consequently, 
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12C(3) ef/^ , 



1 + dt^-'^y 



(38) 



with the recognition that the dimensionless factor e^/"^/(l + e^^'^) is independent of time. 
Since the quantity is of the order of the first term in the sum containing this is 
equivalently an 0{(3^) quantity that is generally negligible. The remainder of the sum may 
be similarly ignored as higher order time derivatives of are successively smaller by factors 
of jS"^, leading to further simplification of Eq. ( ^6|) to 

- f{D-iX)dt' /* j''{D-iX)dt" PDK^ 

D + iX 



f3DK^ 



+ o(/?' 



t'=t 



D2 + A2 



and likewise for its complex conjugate in Eq. (|3^) . 

The last real term in the Eq. (|3^ ) may be shown to reduce to 



(39) 



^ f2Ddt' f' r'2Ddt" P^DK, 



2(3DK^ 



D2 + A2 



+ 



N 

E 



-dt' 



A " io ' [20) Z}2 + A2df(iv+i)^-^^ 



2(3DK^ 



D2 + 



(40) 



t'=t 



by similar arguments such that Eqs. (pGf ) and (|40|) combine to produce 



Pyit) 



Z}2 + A2 



P,(t) + 0{I3'). 



(41) 



Equation (^T]) clearly shows the miraculous cancellation between terms pertaining to a 
nonzero repopulation function originally present in Eq. (|3^ to first order in (5. The fi- 
nal expression for Py{t) is therefore identical to the standard adiabatic result given by Eq. 
(p!3|) in the limit of interest. 

Our next task is to verify that the assumption of instantaneous repopulation in Eq. (plSf ) 
is indeed correct. For this purpose, we begin by expressing as a function of Pq by way of 
Eq. (P, 



(42) 



To the lowest order in /?, Eq. (^2|) is equivalently 



Jo L*^ + A"^ 



(43) 



9 



on which we apply the same technique of integration by parts in conjunction with the 
rationale employed previously to eliminate the various time derivatives to obtain 



P,(t)~-Po(t)+2K,(t) + 2 



■ N 

E 



{\2DJ dfi"") " 



t'=t 



-Po{t)+2K^{t) + 0{P^). 



\2d) D^ + X^dt'(^+^) " 



The last approximate equality in Eq. 



N'^'iip, 0) 



is identically 

A^^q(p, 0)^ 



^ K — 



(44) 



(45) 



which clearly affirms the validity of the instantaneous repopulation assumption to the lowest 
order in f3. 



V. CONCLUSION 



We have taken the exact four-component QKEs for a two-flavour active-sterile neutrino 
system and derived from them an approximate evolution equation for the relic neutrino 
asymmetry in which the role of repopulation is methodically embedded. The consequence 
of including a finite refilling function is to generate higher order terms which are readily 
discarded in the high temperature epoch of interest, thereby yielding a rate equation identical 
to that found earlier in the standard adiabatic limit where the said function is taken to be 
neghgible. The formal procedure developed in the present work to establish this result has 
been labelled the extended adiabatic approximation. 

Frequently adopted in numerical studies with notable accuracy, the assumption of instan- 
taneous repopulation is also shown to arise naturally in the extended adiabatic limit. We 
have thus furnished a rigorous justification for the superficially incompatible assumptions of 
a negligible refilling function and instantaneous repopulation in the regime where collisions 
are the dominant asymmetry generation mode. 
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